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! Abstract 

CN . vrvr scattering length in the 1 = 2 channel is calculated within quenched approxi- 
mation using improved gauge and improved Wilson fermion actions on anisotropic 

^ ' lattices. The results are extrapolated towards the chiral, infinite volume and con- 

■ tinuum limit. This result improves our previous result on the scattering length. In 
. the chiral, infinite volume and continuum limit, we obtain Gq = —0.0467(45), 

Tj" I which is consistent with the result from Chiral Perturbation Theory, the experiment 

■ and results from other lattice calculations. 
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1 Introduction 



Low-energy tttt scattering experiment is a good testing ground for our under- 
standing of the low-energy structure of Quantum Chromodynamics (QCD). 
Chiral Perturbation Theory [1] , Roy equations [2] , dispersion relations [3] and 
other theoretical methods have been used in the study of low-energy tttt scat- 
tering. Lattice QCD is a genuine non-perturbative method which can han- 
dle hadron-hadron scattering at low-energies. In recent years, owing to the 
progress in computing facilities and the use of improved actions, there have 
been several studies on pion pion [4,5,6,7,8,9], pion nucleon [10], pion kaon [11] 
and kaon nucleon [12] scattering in quenched lattice QCD. Very recently, CP- 
PACS collaboration even calculated pion pion scattering phase shift in the 
1 = 2 channel using their unquenched configurations [13,14]. In this letter, we 
report our recent quenched lattice result on the tttt scattering length in the 
1 = 2 channel. The result in this letter is an improvement over our previous 
result on the scattering length [5]. 
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The lattice gauge action used in our study is the tadpole improved action on 
anisotropic lattices [15,16]. Using this gauge action together with improved 
Wilson fermion actions, glueball and hadron spectra have been studied within 
quenched approximation [15,16,17,18,19,20,21]. Configurations generated from 
the gauge action have also been utilized to calculate the the scattering length 
for TTTT scattering in the 1 = 2 channel [5], KN scattering in the 1 = 1 chan- 
nel [12] and the Kn scattering in the 7 = 3/2 channel [11]. In this letter, 
we update our result on the tttt scattering length in the 1 = 2 channel. The 
number of gauge field configurations has been increased. A more general chiral 
extrapolation method is used. Otherwise, the basic procedure of the calcula- 
tion is similar to that adopted in Ref. [5]. 



2 Numerical calculation of the scattering length 

Configurations used in this calculation are generated using the pure gauge 
action for 6^40, 8^40 and 10^50 lattices with the gauge coupling P = 1.9, 2.2, 
2.4, 2.6 and 3.0. The spatial lattice spacing is roughly between O.lfm and 
0.4fm while the spatial physical size of the lattice ranges from 0.7fm to 4.0fm. 
For each set of parameters, several hundred (typically 512) decorrelated gauge 
field configurations are used to measure the fermionic quantities. Statistical 
errors are analyzed using the jack-knife method. Single pion, and rho mass 
values are obtained from the plateau of their corresponding effective mass 
plots with the fitting interval being automatically chosen by minimal per 
degree of freedom. 

In order to calculate the elastic scattering lengths for hadron-hadron scattering 
on the lattice, or the scattering phase shifts in general, one uses Liischer's 
formula which relates the exact energy level of two hadron states in a finite 
box to the elastic scattering phase shift in the continuum. For tttt scattering at 
zero relative three momentum, this formula amounts to a relation between the 
exact energy E^^J of the two pion system with vanishing relative momentum in 
a finite box of size L with isospin /, and the corresponding scattering length 
Qq^ in the continuum. This formula reads [22]: 



where Ci = —2.837297, C2 = 6.375183 are numerical constants and m^^ is the 
mass of the pion. In this letter, we focus on the tttt scattering length a^^ in 
the isospin 1 = 2 channel. 

To measure the hadron mass values ttIt^, rup and to extract the energy shift 
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SE^J, one constructs the correlation functions from the corresponding one 
meson and two meson operators in the appropriate symmetry channel. In this 
letter, we used the same operators as those in Ref . [5] . Numerically, it is more 
advantageous to construct the ratio of the correlation functions: 

n'=\t) = c{,f{t)/(c^(t)c^(t)) , (2) 



where C^^^(t) is the two-pion correlation function in the 1 — 2 channel and 
C7r(t) is the one-pion correlation function. One then uses the fitting function: 

n'=^{t) oc e-'""'^^' - 1 - 6Ei^^t , (3) 



to determine the energy shift SE^^. Usually, the linear form is sufficient. How- 
ever, for small lattices, the exponential form should be used. Similar situation 
was also observed by the JLQCD and CP-PACS collaboration in their calcu- 
lations [6,9]. 

Two pion correlation function, or equivalently, the ratio 71"^^^ (t) is constructed 
from quark propagators, which arc obtained using the Multi-mass Minimal 
Residue algorithm with wall sources. Periodic boundary condition is applied to 
all three spatial directions while in the temporal direction, Dirichlet boundary 
condition is utilized. 

After obtaining the energy shifts SE^^^ from the two-pion correlation functions, 
the values of SE^Jaie substituted into Liischer's formula to solve for the scat- 
tering length Qq'^ for all parameter sets that have been simulated. Prom these 
results, attempts are made to perform an extrapolation towards the chiral, 
infinite volume and continuum limit. 



3 Extrapolations of the scattering length 



The chiral extrapolations of physical quantities are performed as discussed 
in Ref. [5]. The only difference is that, for all physical quantities, since we 
now have more bare quark mass values available, we tried to use a quadratic 
extrapolation in the bare quark mass rather than a linear extrapolation as used 
in Ref. [5] . For each set of parameters, we have calculated fermionic quantities 
with 8 different values of k. The largest value of k, which corresponds to 
lowest pion mass for each simulation point is chosen such that the solution of 
the quark propagator can be obtained within a reasonable number of Minimal 
Residue iterations (typically 600 to 800). The lowest pion to rho mass ratio 
for our simulation points are between 0.66 and 0.74. For the scattering length 
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ao^\ we use the quantity [5] F = ao^^m^/mjr, which in the chiral hmit is given 
by the current algebra value: [23] 



~ -1.364 . (4) 



Here the final numerical value is obtained by substituting in the experimental 
values for rup ~ 770MeV and /jr ~ 93MeV. Chiral perturbation theory with 
loops yields almost the same numerical value for factor F [2] in the 1 = 2 
channel. Note that the factor F can be calculated on the lattice with good 
precision without the lattice calculation of meson decay constants. We have 
adopted a quadratic functional form in chiral extrapolation for all physical 
quantities and the fitting range of the extrapolation is self-adjusted by the 
program to yield a minimal P^r degree of freedom. 
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Fig. 1. Infinite volume extrapolation for tlic quantity F = Qq vnjj/m-,^ obtained from 
our simulation results at /3 = 3.0, 2.6, 2.4, 2.2 and 1.9. The straight line represents 
the linear extrapolation in 1/L^ . The extrapolated results are shown as solid square 
with the corresponding error at = 0. 

The values of F after the chiral extrapolations are used to extrapolated to- 
wards the infinite volume limit by a function linear in 1/L^, as suggested by 
Liischer's formula. This extrapolation is shown in Fig 1 for all five values of 
(3. The extrapolated results thus obtained for the factor F are then used for 
the continuum limit extrapolation. 
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Fig. 2. Continuum extrapolation for the quantity F = a'o'^m'j^/m-^ obtained from our 
simulation results at /? = 3.0, 2.6, 2.4, 2.2 and 1.9. The straight line represents the 
linear extrapolation in a^/ro. The extrapolated result is also shown as solid square 
with the corresponding error. Also shown near Ug = are the (two-loop) result from 
chiral perturbation theory [2] (open hexagon) and the experimental result (shaded 
region). 

Finally, an continuum limit extrapolation is performed to eliminate the finite 
lattice spacing errors. Since we have used the tadpole improved clover Wilson 
action, all physical quantities differ from their continuum counterparts by 
terms that are proportional to the spatial lattice spacing ag- The physical value 
of as in terms of the hadronic scale tq for each value of f3 can be found from Ref. 
[16,21]. The result of the continuum extrapolation is shown in Fig. 2 where 
the results from the chiral and infinite volume extrapolation discussed above 
are indicated as data points (open squares with error bars) in the plot for all 5 
values of /? that have been simulated. The straight line shows the extrapolation 
towards the = limit and the final extrapolated result is also shown as a 
solid square. For comparison, the corresponding result from two-loop chiral 
perturbation theory [2] is also shown near = as the open hexagon. The 
current algebra result and the one-loop chiral perturbation result are quite 
close to the two-loop chiral result numerically. To avoid crowdedness, they are 
not shown in the figure. The experimental result from E865 Collaboration [24] 
is shown as the shaded region near = 0. The height of the shaded region 
designates the error for the experimental result. It is evident that these results 
are compatible within error bars. 

To summarize, we obtain from the continuum extrapolation the following re- 
sult: F = ao'^mp/m^ = —1.41(14). If we substitute in the physical meson 
mass values, we obtain the tttt scattering length in the 1 — 2 channel as: 
af^rriT^ — —0.0467(45), which is to be compared with the current algebra value 
of —0.046 [23], one-loop chiral perturbation theory result of —0.042 [1] and the 
two-loop result of —0.0444(10) [2]. Our result is also consistent with the exper- 
imental result from E865 Collaboration which is: ag ra.,^ = —0.036(9) [24]. The 
result reported in this letter is also consistent with our previous result with 
the error reduced by a factor of two or so due to higher statistics and lattices 
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with finer lattice spacing. This result also agrees with other lattice studies 
reported recently [6,9]. Their final result for Oo^tt is: % tUt^ — —0.0410(69), 
which is also compatible with our result within error. 



4 Conclusions 

In this letter, we report our lattice result on pion-pion scattering lengths in 
isospin 1 = 2 channel, obtained using quenched lattice QCD. Simulations are 
performed on lattices with various sizes, ranging from 0.7fm to about 4fm 
and with five different values of lattice spacing with several hundred gauge 
field configurations. Quark propagators are measured with 8 different valence 
quark mass values. These enable us to explore the chiral limit, the finite volume 
and the finite lattice spacing errors in a systematic fashion. The lattice result 
for the scattering length is extrapolated towards the chiral, infinite volume 
and continuum limit where a result consistent with the experiment, Chiral 
Perturbation Theory and previous lattice results is found. 
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